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Abstract 

The emergence of fractal features in the microscopic structure of space-time is a common theme in 
many approaches to quantum gravity. In this work we carry out a detailed renormalization group 
study of the spectral dimension ds and walk dimension associated with the effective space-times 
of asymptotically safe Quantum Einstein Gravity (QEG). We discover three scaling regimes where 
these generalized dimensions are approximately constant for an extended range of length scales: 
a classical regime where dg = d,dw = 2, a semi-classical regime where dg — 2d/(2> -\- d^,dri^ — 
2 + d, and the UV-fixed point regime where ds = d/2,dw = 4. On the length scales covered 
by three-dimensional Monte Carlo simulations, the resulting spectral dimension is shown to be 
in very good agreement with the data. This comparison also provides a natural explanation 
for the apparent puzzle between the short distance behavior of the spectral dimension reported 
from Causal Dynamical Triangulations (CDT), Euclidean Dynamical Triangulations (EDT), and 
Asymptotic Safety. 



1 Introduction 



Shortly after its introduction, it became clear that the average action approach for gravity [T] 
provides strong support for Weinberg's Asymptotic Safety conjecture [2l|3]. The key ingredient 
in this scenario is a non-trivial fixed point of the gravitational renormalization group flow which 
controls the UV-behavior of the theory and renders Quantum Einstein Gravity (QEG) a non- 
perturbative renormalizable quantum field theory [1,4-33], see |34 1 1351 [36| I37j for reviews. Already 
at a very early stage of the program various indications pointed in the direction that in this 
theory space-time should have certain features in common with a fractal. In ref. [6] the four- 
dimensional graviton propagator has been studied in the regime of asymptotically large momenta 
and it has been found that near the Planck scale a kind of dynamical dimensional reduction 
occurs. As a consequence of the non-Gaussian fixed point (NGFP) controlling the UV behavior 
of the theory, the four-dimensional graviton propagator essentially behaves two-dimensional on 
microscopic scales. 

Subsequently, the finger prints of the NGFP on the fabric of the effective QEG space-times 
have been discussed in [TU], where it was shown that asymptotic safety induces a characteristic self- 
similarity of space-time on length-scales below the Planck length £pl. The graviton propagator 
becomes scale- invariant in this regime 0, and based on this observation it was argued that, 
in a cosmological context, the geometry fluctuations it describes can give rise to a scale free 
spectrum of primordial density perturbations responsible for structure formation |38t 139) . Thus 
the overall picture of the space-time structure in asymptotically safe gravity as it emerged about 
ten years ago comprises a smooth classical manifold on large distance scales, while on small 
scales one encounters a low dimensional effective fractal [6l[T0]. The feature at the heart of these 
results is the observation that the effective field equations derived from the gravitational average 
action equip every given smooth space-time manifold with, in principle, infinitely many different 
(pseudo) Riemannian structures, one for each coarse graining scale [HI [15]. Thus, very much like 
in the famous example of the coast line of England [40] , the proper length on a QEG space-time 
depends on the length of the yardstick used to measure it0 

Along a different line of investigations, the Causal Dynamical Triangulation (CDT) approach 
has been developed and first Monte-Carlo simulations were performed |43j . see [13| for a recent 
review. In this framework one attempts to compute quantum gravity partition functions by 
numerically constructing the continuum limit of an appropriate statistical mechanics system. 
From the perspective of the latter, this limit amounts to a second order phase transition. If CDT 
and its counterpart QEG, formulated in the continuum by means of the average action, belong to 
the same universality clasj§ one may expect that the phase transition of the former is described 
by the non-trivial fixed point underlying the asymptotic safety of the latter. 

Remarkably, ref. [l5] reported results which indicated that the four-dimensional CDT space- 
times, too, undergo a dimensional reduction from four to two dimensions as one "zooms" in on 

^Earlier on similar fractal properties had already been found in other quantum gravity theories, in particular 
near dimension 2 [41] and in a non-asymptotically safe model [42| . 

^For the time being this is merely a conjecture, of course, albeit a very natural one. 
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short distances. In particular it had been demonstrated that the spectral dimension dg measured 
in the CDT simulations has the very same limiting behaviors, 4 — t- 2, as in QEG [46J. Therefore 
it was plausible to assume that both approaches indeed "see" the same continuum physics. 

However, this interpretation became problematic when it turned out that the Monte Carlo 
data correspond to a regime where the cutoff length inherent in the triangulations is still signifi- 
cantly larger than the Planck length. According to the renormalization group (RG) trajectories 
computed in QEG one would not expect that the asymptotic scaling behavior implied by the fixed 
point is already realized there [9]; on the other hand it is exactly this asymptotic scaling regime to 
which the QEG prediction of 2 fractal dimensions pertains [6l \TU[ H6] . Thus the obvious question 
is why the CDT simulations detect a significant dimensional reduction despite their appreciable 
distance to the continuum limit. 

Recently the situation became even more puzzling. In particular, ref. |47] carried out CDT 
simulations for d = 3 macroscopic dimensions, which favor a value near = 2 on the shortest 
length-scale probed; in this case the QEG prediction for the fixed point region is the value 
ds = 3/2, however Furthermore, the authors of ref. [35] reported simulations within the 

euclidean dynamical triangulation (EDT) approach in (i = 4, which favor a drop of the spectral 
dimension from 4 to about 1.5; this is again in conflict with the QEG expectations if one interprets 
the latter dimension as the value in the continuum limit. 

One of the aims of the present paper is to propose a resolution to these puzzles. In this course 
we will explicitly compute several types of scale dependent effective dimensions, specifically the 
spectral dimension ds and the walk dimension d^ for the effective QEG space-times. We shall 
see that on length scales slightly larger than there exists a further regime which exhibits the 
phenomenon of dynamical dimensional reduction. There the spectral dimension is even smaller 
than near the fixed point, namely = 4/3 in the case of 4 dimensions classically. Moreover, we 
shall demonstrate in detail that the (3-dimensional) results reported in |47j are in perfect accord 
with QEG. In this course, we also verify the supposition [17j that the shortest possible length 
scale achieved in the simulations is not yet close to the Planck length. Rather the Monte Carlo 
data probes the transition between the classical and the newly discovered "semi-classical" regime. 

It is intriguing that Loop Quantum Gravity and spin foam models also show indications for 
a similar dimensional reduction |49^ [50] , with some hints for an intermittent regime where the 
spectral dimension is smaller than in the deep ultraviolet. In ref. |51) an argument based upon the 
strong coupling limit of the Wheeler- De Witt equation was put forward as a possible explanation 
of this dimensional reduction. Within non-commutative geometry Connes et al. [52] interpreted 
the dynamical dimensional reduction lo dg = 2 which was observed in QEG in the context of 
the derivation of the Standard Model from a spectral triple. In fact, from the data provided by 
a spectral triple, its Dirac operator in particular, one can compute a type of spectral dimension 
of the resulting non-commutative space which is closely related to the one we are considering 
here. Also for standard fractals such as Cantor sets, it has been possible to find spectral triples 
representing them and to compute the corresponding dimensions [53J. 

Furthermore, a number of model systems (quantum sphere, K-Minkowski space, etc.) give 
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rise to a similar reduction as quantum gravity |54j . Among other developments, these findings 
also motivated the investigation of physics on prescribed fractal space-times. In refs. [55\ [56] 
a fractional differential calculus [57j was employed in order to incorporate fractal features, and 
in [58] recent exact results on spectral zeta-functions on certain fractals [59] were used to study 
the thermodynamics of photons on fractals. In ref. [60] matter quantum field theories were 
constructed and renormalized on a fractal background. This almost universal appearance of 
fractional properties of space-time and its accessibility in various, a priori different, approaches 
to Quantum Gravity make the generalized notions of dimensionality discussed in this paper a 
valuable tool in comparing the physics content of these different formulations. 

The remaining parts of this paper are organized as follows. In Sect. [2] we introduce the 
different notions of "dimension", i.e., the spectral, walk, and Hausdorff dimension which will be 
used to characterize the fractal properties of the effective QEG space-times in Sect. [3l In Sect. H] 
we then analyze how these dimensions change with the RG-scale and identify the scaling regimes 
where they are approximately constant. These results are then compared with the CDT data 
obtained in [17] in Sect. [SJ We conclude with a brief discussion of our findings in Sect. 

2 Generalized dimensions characterizing fractal space-times 

Investigating random walks and diffusion processes on fractals, one is led to introduce various 
notions of fractal dimensions, such as the spectral or walk dimension [61]. These notions also 
prove useful when characterizing properties of space-time in quantum gravity, and we will review 
these concepts in the remainder of this section. 

2.1 The spectral dimension 

To start with, consider the diffusion process where a spin-less test particle performs a Brownian 
random walk on an ordinary Riemannian manifold with a fixed classical metric g^i,{x). It is 
described by the heat-kernel Kg{x,x';T) which gives the probability density for a transition of 
the particle from x to x' during the fictitious time T. It satisfies the heat equation 

drKgix, x'; T) = -AgKgix, x'; T) , (2.1) 

where = —D^ denotes the Laplace-Beltrami operator. In flat space, this equation is easily 
solved by 

Kg{x,x'-T) = I ^e^^ (— ')e-^^^ (2.2) 

In general, the heat-kernel is a matrix element of the operator exp(— TA^). In the random walk 
picture its trace per unit volume, 

Pg{T) = J d'^x./^Kgix, X- T) = Tr exp(-rAg) , (2.3) 

has the interpretation of an average return probability. Here V = j d'^x^J g{x) denotes the total 
volume. It is well known that Pg possesses an asymptotic early time expansion (for T — )• 0) of 
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the form Pg{T) = (AnTf^'^ E^o^n^"> with denoting the Seeley-DeWitt coefficients. From 
this expansion one can motivate the definition of the spectral dimension ds as the T-independent 
logarithmic derivative 

d 



d\nPg{T) 



dlnT 



(2.4) 

T=0 

On smooth manifolds, where the early time exapnsion of Pg{T) is valid, the spectral dimension 
agrees with the topological dimension d of the manifold. 

Given Pg(T), it is natural to define an, in general T-dependent, generalization of the spectral 
dimension by 

_ dlnP„(T) , , 

According to (|2.4p . we recover the true spectral dimension of the space-time by considering the 
shortest possible random walks, i.e., by taking the limit ds = limr_>.o (T) . Note that in view 
of a possible comparison with other (discrete) approaches to quantum gravity the generalized, 
scale-dependent version (12. 5p will play a central role later on. 



2.2 The walk dimension 

Regular Brownian motion in fiat space has the celebrated property that the random walker's 
average square displacement increases linearly with time: (r^) cx T. Indeed, performing the 
integral ()2.2p we obtain the familiar probability density 

K{x,x'-T) = {4nT)-'/'e^p {-^^) (^.6) 

with a{x, x') = i|x — x'p half the squared geodesic distance between the points x, x' . Using (|2.6p 
yields the expectation value (r^) = (x^) = / d'^x x"^ K{x,^]T) oc T. 

Many diffusion processes of physical interest (such as diffusion on fractals) are anomalous in 
the sense that this linear relationship is generalized to a power law (r^) oc T'^I'^^ with dw ^ '2,. 
The interpretation of the so-called walk dimension dw is as follows. The trail left by the random 
walker is a random object, which is interesting in its own right. It has the properties of a fractal, 
even in the "classical" case when the walk takes place on a regular manifold. The quantity dr^ is 
precisely the fractal dimension of this trail. Diffusion processes are called regular if dw = 2, and 
anomalous when dyj ^2. 



2.3 The Hausdorff dimension 

Finally, we introduce the Hausdorff dimension dn- Instead of working with its mathematical 
rigorous definition in terms of the Hausdorff measure and all possible covers of the metric space 
under consideration, the present, simplified definition may suffice for our present purposes. On a 
smooth set, the scaling law for the volume V{r) of a d-dimensional ball of radius r takes the form 

V{r) oc r'^" . (2.7) 
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The Hausdorff dimension is tlien obtained in the limit of infinitely small radius, 

InV(r) , , 

du = lim . ^ ' . (2.8 
r-i-o Inr 

Contrary to the spectral or walk dimension whose definitions are linked to dynamical diffusion 
precesses on space-time, there is no dynamics associated with du- 

3 Fractal dimensions within QEG 

Upon introducing various concepts for fractal dimensions in the last section, we now proceed 
with their evaluation for the QEG effective space-times, generalizing the results of ref. [IB]. Our 
discussion will mostly be based on the so-called Einstein-Hilbert truncation introduced in the 
next subsection. As we shall see this restriction is actually unnecessary in the asymptotic scaling 
regime, i.e., when the RG-trajectory is close to the NGFP. In this case we can derive exact results 
for the spectral and walk dimension by exploiting the scale invariance of the theory at the fixed 
point. 

3.1 Diffusion processes on QEG space-times 

Since in QEG one integrates over all metrics, the central idea is to replace Pg{T) by its expectation 
value 

P{T) = {P^{T)) = J p7PCPC'P^(r)e~^ba.c[7,c,c] ^ ^3^;^) 

Here 7^,^ denotes the microscopic metric and 5bare is the bare action related to the UV fixed 
point, with the gauge- fixing and the pieces containing the ghosts C and C included. For the 
untraced heat-kernel, we define likewise 

K{x,x';T) = {K^{x,x';T)) . (3.2) 

These expectation values are most conveniently calculated from the effective average action Tk, 
which equips the d-dimensional smooth manifolds underlying the QEG effective space-times with 
a family of metric structures {{g^iv)k-,^ < < 00}, one for each coarse- graining scale k [lUJIH]. 
These metrics are solutions to the effective field equations implied by r^. 

To start with, we shall approximate the latter by the Einstein-Hilbert truncation [H [9] 

r. = (16,0.)-' / <i^.v^ i-R + 2A.) + classical gauge-Axing and ghost term. , (3.3) 

which includes a scale-dependent cosmological constant and Newtons constant Gk- The cor- 
responding effective field equation reads 

2 - 

It has the same form as the classical Einstein equation, with a /c-dependent cosmological constant 
Afc, however. We can easily find the fc-dependence of the corresponding solution {g^u)k by rewrit- 
ing (|3.4p as [Xko/^k]R'^u{{g)k) = ^r^^ko^^u for some fixed reference scale ko, and exploiting that 
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R^v[cg) = c ^R^y{g) for any constant c > 0. This shows that the metric and its inverse scale 
according to 

{9A^))k = [KJM{9A^))k, , {g^''{x))k = [Xk/~XkoW{x))k„ ■ (3.5) 

Denoting the Laplace-Beltrami operators corresponding to the metrics {g^u)k and {g^u)ko by A(A;) 
and A(/co)) respectively, these relations imply 

A(A;) = [Xk/K] A(A;o) . (3.6) 

At this stage, the following remarks are in order. In the asymptotic scaling regime associated 
with the NGFP the scale-dependence of the couplings is fixed by the fixed point condition: 

Afc oc A;^ Gfc oc k^'^ . (3.7) 

This implies in particular 

{gf^uix))k oc k""^ (k^oo). (3.8) 

This asymptotic relation is actually an exact consequence of Asymptotic Safety, which solely relies 
on the scale-independence of the theory at the fixed point. 

In general, the relation p.6p will receive corrections from the presence of higher-derivative 
operators in the effective average action [10 \ \16\ [T8\ I22|. These contributions organize themselves 
into a power series controlled by the dimensionless cosmological constant Xk = X'l/k'^ in which 
the square bracket in (|3.6p provide the leading term for small values of Afc. Thus, for finite scales 
k, the equations (j3.5p and (j3.6p hold true in the Einstein-Hilbert truncation and the absence of 
matter only. These conditions guarantee that the effective field equation has the simple form 
(j3.4p . which is necessary for the proportionality {gfj,y)k oc A^^ expressed in (j3.7p . This relation 
constitutes an essential piece in the derivation of eq. (|3.6p . In the scaling regime of a complete 
gravity-matter fixed point the conditions above can be relaxed, however, since then for purely 
dimensional reasons basically, {g^u)k oc in the full theory and any sensible truncation. 

This said, we can now evaluate the expectation value (j3.ip by exploiting the effective field 
theory properties of the effective average action. Since defines an effective field theory at 
the scale k we know that {0{^^y)) 0{{g^u)k) provided the observable O involves only mo- 
mentum scales of the order of k. We apply this rule to the RHS of the diffusion equation, 
O = —A^K^{x,x';T). The subtle issue here is the correct identification of k. If the diffusion 
process involves (approximately) only a small interval of scales near k over which A^ does not 
change much, the corresponding heat equation contains the operator A(/c) for this specific, fixed 
value of k: dTK{x,x';T) = —A{k)K{x,x';T). Denoting the eigenvalues of A(A;o) by and the 
corresponding eigenfunctions by (pn, this equation is solved by 

4>n{x)(t>n{x') exp - F{k'^)£nT^ ■ (3.9) 

n 

Here we introduced the convenient notation F{k'^) = Afc/Afc,,. Knowing the propagation kernel, 
we can time-evolve any initial probability distribution p{x; 0) according to 

p{x;T)= [ d''xWgoix')K{x,x';T)p{x';0) (3.10) 
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with go the determinant of {gp,y)ka- If the initial distribution has an eigenfunction expansion of 
the form p{x; 0) = ^„ Cn(t>n{x) we obtain 

p{x- T) = Y, CnMx) exp ( - F(A;2)£:„r) . (3.11) 

n 

If the CnS are significantly different from zero only for a single eigenvalue Sn, we are dealing 
with a single-scale problem and would identify k'^ = £]\f as the relevant scale at which the running 
couplings are to be evaluated. In general the C^'s are different from zero over a wide range 
of eigenvalues. In this case we face a multiscale problem where different modes (pn probe the 
space-time on different length scales. If A(A;o) corresponds to flat space, say, the eigenfunctions 
4'n = (t>p are plane waves with momentum p'^, and they resolve structures on a length scale ^ of 
order l/|p|. Hence, in terms of the eigenvalue £n = £p = the resolution is £ l/\f£^. This 
suggests that when the manifold is probed by a mode with eigenvalue it "sees" the metric 
{g^u)k for the scale k = \fE^. Actually, the identification k = \fE^ is correct also for curved 
space since, in the construction of F^, the parameter k is introduced precisely as a cutoff in the 
spectrum of the covariant Laplacian. 

As a consequence, under the spectral sum of (13. lip , we must use the scale fc^ = En which 
depends explicitly on the resolving power of the corresponding mode. Likewise, in eq. (j3.9p . F{k'^) 
is to be interpreted as F{£n)- 

K{x, x'; T) = ^ (pn{x)(j)n{x') exp - F{£n)£nT^ 

" (3.12) 
= ^(/)„(x')exp(^-F(A(A;o))A(A:o)r),/. 

n (^ ) • 

n 

As in [1^, we choose k^ as a macroscopic scale in the classical regime, and we assume that at /cq 
the cosmological constant is small, so that {g^v)k(, can be approximated by the flat metric on W^. 
The eigenfunctions of A(/co) are plane waves then and eq. (j3.12p becomes 

K{x, x'- T)= [ e'P<--^'^ e-^'np')T (3.13) 

J {^T^r 

where the scalar products are performed with respect to the flat metric, {gnu)ko = ^f^u- The 
kernel ()3.13p satisfies K{x,x';0) = 6'^{x — x') and, provided that \\mp^Qp^F{p'^) = 0, also 
/ d'^xK{x,x';T) = 1. 

Taking the normalized trace of ()3.13p within this "fiat space-approximation" yields [l6] 
Introducing z = p^, the final result for the average return probability reads 

feo- 



where F{z) = \{k'^ = z)/\ 
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3.2 The spectral dimension in QEG 

In the classical case, F{z) = 1, the relation (13.15P reproduces the familiar result P{T) = 
l/(47^T)'^/^ whence Vs{T) = d independently of T. We shall now discuss the spectral dimension 
for several other illustrative and important examples. 

(A) To start with, let us evaluate the average return probability (13.151) for a simplified RG- 
trajectory where the scale dependence of the cosmological constant is given by a power law, with 
the same exponent 5 for all values of k: 

Afc oc fc^ =^ F{z) oc (3 J6) 

By rescaling the integration variable in (jS.lSp we see that in this case 



nT) = ^. (3.17) 



Hence (|2.5p yields the important result 



(3.18) 



It happens to be T-independent, so that for T — )• trivially 



4 = 5^. (3.19) 

(B) Next, let us be slightly more general and assume that the power law (I3.16P is valid only 
for squared momenta in a certain interval, G [zi,Z2], but remains unspecified otherwise. 
In this case we can obtain only partial information about P{T), namely for T in the interval 
[^2^"*^, The reason is that for T G [z^"*^, Zj""*^] the integral in (|3.15p is dominated by momenta 
for which approximately l/p^ T, i.e., z G [2:1,^2] • This leads us again to the formula (I3.18p . 
which now, however, is valid only for a restricted range of diffusion times T; in particular the 
spectral dimension of interest may not be given by extrapolating (|3.18p to T — )■ 0. 

(C) Let us consider an arbitrary asymptotically safe RG-trajectory so that its behavior for A; — t- 00 
is controlled by the NGFP. In this case the running of the cosmological constant for k > M, with 
M a characteristic mass scale of the order of the Planck mass, is given by a quadratic scale- 
dependence Afc = A*/c^, independently of d. This corresponds to a power law with 6 = 2, which 
entails in the NGFP regime, i.e., for T < 

VsiT) = ^ (nGFF regime) . (3.20) 

This dimension, again, is locally T-independent. It coincides with the T — )• limit: 

ds = ^. (3.21) 

This is the result first derived in ref. [l6]. As it was explained there, it is actually an exact 
consequence of Asymptotic Safety which relies solely on the existence of the NGFP and does not 
depend on the Einstein-Hilbert truncation. 
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(D) Returning to the Einstein-Hilbert truncation, let us consider the piece of the Type Ilia 
RG-trajectory depicted in Fig. [1] which hes inside the hnear regime of the Gaussian fixed point. 
Newton's constant is approximately /c-independent there and the cosmological constant evolves 
according to 

Afc = Ao + lyGok'^. (3.22) 

Here = (47r)^^'^/^((i — 3)<I>^^2(0) is a scheme-dependent constant [BIS]. When k is not too small, 
so that Ao can be neglected relative to vG^k'^, we are in what we shall call the "A;'^ regime"; it is 
characterized by a pure power law A^ ~ k^ with 6 = d. The physics behind this scale dependence 
is simple and well-known: It represents exactly the vacuum energy density obtained by summing 
up the zero point energies of all field modes integrated out. For T in the range of scales pertaining 
to the k'^ regime we find 

= ^ (^'^ ^^Sime) . (3.23) 

Note that for every d > 2 the spectral dimension in the k'^ regime is even smaller than in the 
NGFP regime 

^^(NGFP regime) 



T)s{k'^ regime) 
3.3 The walk dimension in QEG 



l + (d-2)/4. (3.24) 



In order to determine the walk dimension for the diffusion on the effective QEG space-times we 
return to eq. (j3.13p for the untraced heat-kernel. We restrict ourselves to a regime with a power 
law running of A^, whence F{p'^) = (Lp)^ with some constant length-scale L. 

Introducing = p^T^^^'^~^^^ and .^^ = (x^ — x'^ /T^l^'^^^^ we can rewrite (j3.13p in the form 

^(^'^''^) = ^^^(|^) (3-25) 

with the function 

For (5 = 0, this obviously reproduces (|2.6p . From the argument of <I> in 
r = \x — x'\ scales as T^/'y^^^') so that the walk dimension can be read off aa 



3.25P we infer that 



V^{T) = 2 + 5 



(3.27) 



In analogy with the spectral dimension, we use the notation T>w{T) rather than to indicate 
that it might refer to an approximate scaling law which is valid for a finite range of scales only. 
For 5 = 0,2, and d we find in particular, for any topological dimension d, 



2 classical regime 

4 NGFP regime (3.^ 
2 + d k'^ regime 



'Cf. eq. (5.18) in ref. gl]. 
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Regimes with all three walk dimensions of (I3.28P can be realized along a single RG-trajectory. 
Notably, the result for the NGFP regime, = 4, is exact in the sense, that it does not rely on 
the Einstein-Hilbert truncation. 

3.4 The HausdorfF dimension in QEG 

The smooth manifold underlying QEG has per se no fractal properties whatsoever. In particular, 
the volume of a d-ball B"^ covering a patch of the smooth manifold of QEG space-time scales as 

ViB'')= [ d'^x^ocirk)'. (3.29) 

Thus, by comparing to eq. ()2.7p . we read off that the Hausdorff dimension is strictly equal to the 
topological one: 



dn = d 



(3.30) 



We emphasize that the effective QEG space-times should not be visualized as a kind of sponge. 
Their fractal-like properties have no simple geometric interpretation; they are not due to a "re- 
moving" of space-time points. Rather they are of an entirely dynamical nature, reflecting certain 
properties of the quantum states the system "space-time metric" can be in. 

For standard fractals the quantities ds, d^, and dn are not independent but are related by 

m 

By combining eqs. (j3.18p . (j3.27p . and (j3.30p we see that the same relation holds true for the 
effective QEG space-times, at least within the Einstein-Hilbert approximation and when the 
underlying RG-trajectory is in a regime with power-law scaling of Afc. For every value of the 
exponent 5 we have 

= f- . 3.32 

The results dn = d, = 2 + 6 imply that, as soon as 6 > d — 2, we have > dn and the 
random walk is recurrent then [61]. Classically {6 = 0) this condition is met only in low dimensions 
d < 2, but in the case of the QEG space-times it is always satisfied in the k'^ regime {6 = d), 
for example. So also from this perspective the QEG space-times, due to the specific quantum 
gravitational dynamics to which they owe their existence, appear to have a dimensionality smaller 
than their topological one. 

It is particularly intriguing that, in the NGFP regime, = 4 independently of d. Hence the 
walk is recurrent {Dw > dn) for d < 4, non-recurrent for d > 4, and the marginal case Dw = dn 
is realized if and only if d = 4, making d = 4 a distinguished value. Notably, there is another 
feature of the QEG space-times which singles out d = 4: It is the only dimensionality for which 
I?s(NGFP regime) = d/2 coincides with the effective dimension dgg = d + rj^^ = 2 derived from the 
graviton propagator [6l Il6] . 
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The relation (|3.32p also has an important implication for a possible relation between the QEG 
effective space-times and those of the CDT approach. The latter have a non-classical Hausdorff 
dimension dn ^ d on microscopic scales, while dn = d in QEG. Hence, by ()3.32p . we cannot 
expect that both Dg and D^i agree between CDT and QEG. If it should turn out that actually 
in some non-classical regime, then and vT"^ are necessarily different there. 



4 The RG-flow of and 

We now proceed by discussing the scale-dependence of the spectral and walk dimension. For this 
purpose, we consider an arbitrary RG-trajectory k i— )• ((7fc,Afc), where = Gkk'^~'^ and = 
Afc/j"^ are the dimensionless counterparts of Newton's constant and the cosmological constant, 
respectively. Along such a RG-trajectory there might be isolated intervals of A;- values where the 
cosmological constant evolves according to a power law, A^ oc k^, for some constant exponents 
6 which are not necessarily the same on different such intervals. If the intervals are sufficiently 
long, it is meaningful to ascribe a spectral and walk dimension to them since 5 = const implies 
/c-independent values = 2d/ (2 + 6) and P^, = 2 + 5. 

In between the intervals of approximately constant and Dw, where the /c-dependence of 
Afc is not a power law, the notion of a spectral or walk dimension might not be meaningful. The 
concept of a scale- dependent dimension or is to some extent arbitrary with respect to the 
way it interpolates between the "plateaus" on which 6 = const for some extended period of RG 
time. While RG methods allow the computation of the T>s and 2?^ values on the various plateaus, 
it is a matter of convention how to combine them into continuous functions k i— >■ 'Ds{k),'Dw{k) 
which interpolate between the respective values. 

4.1 The exponent 5 as a function on theory space 

In this subsection, we describe a special proposal for a /c-dependent T>s{k) and T)y^{k) which is 
motivated by technical simplicity and the general insights it allows. We retain eqs. (j3.18p and 
(j3.27p . but promote 5 ^ 5{k) to a A;-dependent quantity 

5{k) = kdk\n{\k). (4.1) 

When Afc satisfies a power law, A^ oc k^ this relation reduces to the case of constant 5. If not, 

5 has its own scale dependence, but no direct physical interpretation should be attributed to 
it. The particular definition (j4.ip has the special property that it actually can be evaluated 
without first solving for the RG-trajectory. The function 5{k) can be seen as arising from a 
certain scalar function on theory space, 5 = 6{g,X), whose A;-dependence results from inserting 
an RG-trajectory: 5{k) = S{gk,Xk)- In fact, (j4.ip implies 5{k) = kdkln{k'^Xk) = 2 + X^^kd^Xk 
so that 6{k) = 2 + X^^ f3x{gk, X^) upon using the RG-equation kd^Xk = l3x{g,X). Thus when we 
consider 6 as a function on theory space, coordinatized by g and A, it reads 

6{g,X) = 2 + \pxig,X). (4.2) 
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Substituting this relation into ()3.18p and (13.271) . the spectral and the walk dimensions become 
functions on the (7-A-plane 



2d 



4 + A-1/3a(5,A) ' 



and 



P^(5,A) = 4 + A-Va(9,A). 
To evaluate these expressions further, we use the /3-functions derived in |1]: 

/?A(5,A)=(r?;v-2)A + i(47r)i-^/2g 

X 2d{d + l)$i/2(-2A) - 8d$^/2(0) - d{d + l)r?^l>i/2(-2A) 
f3g{g,X) = {d-2 + r]N)g. 
Here the anomalous dimension of Newton's constant r/7v is given by 

1 - 5^2 (A) 

with the following functions of the dimensionless cosmological constant: 

B,{X) = i (4vr)i-'^/2 ^ 1)^1 (_2A) - 6d{d - 1)<^1,(,-2X) 



d/2- 



-i(0)-24<I>^/2(0) 



i?2(A)= -i(47r: 



l-d/2 



-2X)-6d{d-l)^l/,{-2X) 



(4.3) 
(4.4) 



(4.5) 



(4.6) 



(4.7) 



For practical computations we use the threshold functions resulting from the optimized cutoff 

1 1 = 1 1 



KM 



r(n + 1) {l+w)P 



$P(u;) 



r(n + 2) {l+w)P ' 



(4.8) 



As we discussed already, the scaling regime of a NGFP has the exponent 5 = 2. From eq. 
(j4.2p we learn that this value is realized at all points {g, A) where /3a = 0. The second condition 
for the NGFP, f3g = 0, is not required here, so that we have 6 = 2 along the entire line in theory 
space: 

13 = [{g,X)\Px{9,X) = 0}. (4.9) 

For d = 4 the curve B is shown as the bold blue line in Fig. [TJ Both the Gaussian fixed point 
(GFP) (5, A) = (0,0) and the NGFP, (5, A) = (g*,A=K), are located on this curve|l| Furthermore, 
the turning points T of all Type Ilia trajectories are also situated on B, and the same holds for all 
the higher order turning points which occur when the trajectory spirals around the NGFP. This 
observation leads us to an important conclusion: The values 6 = 2 <^=^> Dg = d/2,'Dw = 4 which 
(without involving any truncation) are found in the NGFP regime, actually also apply to all 
points {g, A) G B, provided the Einstein-Hilbert truncation is reliable and no matter is included. 

"'At the GFP Vb = d, = 2, however, since at this point both A = 0, /3a = so that A~^/3a|gfp = —2. 
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Figure 1: The ^-A— theory space with the hne of turning points, i3, and a typical trajectory 
of Type Ilia. The arrows point in the direction of decreasing k. The big black dot indicates 
the NGFP while the smaller dots represent points at which the RG-trajectory switches from 
increasing to decreasing A or vice versa. The point T is the lowest turning point, and C is a 
typical point within the classical regime. For A > 0.4, the RG-flow leaves the classical regime and 
is no longer reliably captured by the Einstein-Hilbert truncation. 

4.2 The spectral and walk dimensions along a RG-trajectory 

We proceed by investigating how the spectral and walk dimension of the effective QEG space-times 
changes along a given RG-trajectory. As discussed above, our interest is in scaling regimes where 
Vs and T>w remain (approximately) constant for a long interval of A;-values. For the remainder of 
this subsection, we will restrict ourselves to the case d = A for concreteness. 
We start by numerically solving the coupled differential equations 



with the /3-functions (14. Sp for a series of initial conditions keeping Ainit = A(A:o) = 0.2 fixed and 
successively lowering g-^-a\t = giko)- The result is a family of RG-trajectories where the classical 
regime becomes more and more pronounced. Subsequently, these solutions are substituted into 
and which give 5init) '^init) and I'io(t; <7init; -^init) in dependence of the RG-time 

t = ln{k) and the RG-trajectory. One can verify explicitly, that substituting the RG-trajectory 
into the return probability (j3.15p and computing the spectral dimension from ()2.4p by carrying 
out the resulting integrals numerically gives rise to the same picture. 

Fig. [3] then shows the resulting spectral dimension, the walk dimension, and the localization 
of the plateau-regimes on the RG-trajectory in the top-left, top-right and lower diagram, respec- 
tively. In the top diagrams, ginit decreases by one order of magnitude for each shown trajectory. 



kdkg{k)=f3g{g{k),Kk)) 



kdkX{k)=p^{g{k),X{k)) 



(4.10) 
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Figure 2: The line of turning points B divides the (7-A-plane in two parts on which Dg and 
are everywhere either larger or smaller than at the NGFP. 



starting with the highest value to the very left. As a central result, Fig. [3] establishes that the 
RG-flow gives rise to three plateaus where I?s(t) and 'Dw{t) are approximately constant: 

(i) For small values k, below t ~ 1.8, say, one finds a classical plateau where Dg = 4,I'«, = 2 for a 
long range of fc- values. Here 6 = 0, indicating that the cosmological constant is indeed constant. 

(ii) Following the RG-flow towards the UV (larger values of t) one next encounters the semi- 
classical plateau where P<j = 4/3, = 6. In this case 6{k) = 4 so that oc on the 
corresponding part of the RG-trajectory. 

(iii) Finally, the NGFP plateau is characterized by Ds = 2, = 4, which results from the scale- 
dependence of the cosmological constant at the NGFP oc k'^ <;=^ 6 = 2. 



At this stage, it is worthwhile to see which parts of a typical RG-trajectory realize the scaling 
relations underlying the plateau- values of Vg and V^. This is depicted in the lower diagram of 
Fig. [3] where we singled out the third solution with gi^n = 10~^ for illustrative purposes. In this 
case the classical plateau is bounded by the points Pi and P2 and appears well below the turning 
point T, while the semi-classical plateau is situated between the points P3 and P4 well above 
the turning point. The NGFP plateau is realized by the piece of the RG-trajectory between P5 
and the NGFP. The turning point T is not situated in any scaling region but appears along the 
crossover from the classical to the semi-classical regime of the QEG space-times. For t < 0, the 
spectral dimension (walk dimension) increases (decreases) rapidly. In this region, the underlying 
RG-trajectory is evaluated outside the classical regime at points A > 0.35. In this region of the 
theory space, the Einstein-Hilbert truncation is no longer trustworthy, so that this rapid increase 
of Dg is most likely an artefact, arising from the use of an insufficient truncation. 
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Figure 3: The t = ln(/i;)-dependent spectral dimension (upper left) and walk dimension (upper 
right) along illustrative solutions of the RG-equations (j4.10p in d = 4. The trajectories develop 
three plateaus: the classical plateau with Dg = ^,T^w = 2, the semi-classical plateau where 
T>s = 4/3, = 6 and the NGFP plateau with Dg = 2,P^ = 4. These plateau values are 
indicated by the gray horizontal lines and connected by crossover parts. The lower figure shows 
the location of these plateaus on the RG-trajectory: the classical, /c^, and NGFP regime appear 
between the points Pi and P2, P3 and P4, and above P5, respectively. 

Notably, the plateaus observed above become more and more extended the closer the tra- 
jectories turning point T gets to the GFP, i.e., the smaller the IR value of the cosmological 
constant. The first RG-trajectory with the largest value ginit = 0.1 does not even develop a 
classical and semi-classical plateau, so that a certain level of fine-tuning of the initial conditions 
is required in order to make these structures visible. Interestingly enough, when one matches the 
observed data against the RG-trajectories of the Einstein-Hilbert truncation [631 EH] one finds 
that the "RG-trajectory realized by Nature" displays a very extreme fine-tuning of this sort. The 
coordinates of the turning point are approximately qt ^ \t ^ 10"^*^ and it is passed at the 
scale Ut ~ 10~^*^mpi ~ 10~^eV ~ (10~^mm)~^, so that there will be very pronounced plateau 
structures in this case. 
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Figure 4: Spectral dimension I?^^^(r) determined from random walks on a CDT space-time 
built from N = 200k simplices [37]. 



5 Matching the spectral dimensions of QEG and CDT 

The key advantage of the spectral dimension T>s(T) is that it may be defined and computed within 
various a priori unrelated approaches to quantum gravity. In particular, it is easily accessible in 
Monte Carlo simulations of the Causal Dynamical Triangulations (CDT) approach in d = 4 j35] 
and d = 3 [47J as well as in Euclidean Dynamical Triangulations (EDT) [IS]. This feature allows 
a direct comparison between T)^^^{T) and T)f^^(T) obtained within the discrete approaches 
and 2?^^'^(T) capturing the fractal properties of the QEG effective space-times. In this section 
we will carry out this comparison for d = 3. In particular we shall determine the specific RG- 
trajectory of QEG which, we believe, underlies the numerical data obtained in [U]. In principle, 
it is straightforward to do the same comparison in d = 4. This, however, will require access to 
the detailed Monte Carlo data produced by the four-dimensional CDT or EDT simulations If] 

Let us start by looking into the typical features of the spectral dimension V^^^ (T) obtained 
from the simulations. A prototypical data set showing V^^'^iT) as function of the length of the 
random walk T is given in Fig. HI The resulting curve is conveniently split into three regimes: 

(i) For T < 20, corresponding to the left gray region in Fig. U ^^^^{T) undergoes rapid oscilla- 
tions. These originate from the discrete structure of the triangulation to which the short random 
walks are particular sensitive. 

(ii) For long random walks with T > 500, the data shows an exponential fall-off. This feature 
is due to the compact nature of the triangulation, which implies that for long random walks 
pCDT^2^^ is governed by the lowest eigenvalue of the Laplacian on the compact space. This 
regime is marked by the right gray region in Fig. HI 

(iii) Between these two regimes, Pp^^(T) is affected neither by the discreteness nor the com- 
pactness of the triangulation. Since for pP^'~^(T), determined by the flat-space approximation 
p.l4p . we do not expect any of these effects to appear, we use this middle region to compare the 
T-dependent spectral dimensions arising from the two, a priori different, approaches. 

This comparison is then carried out as follows: 

^We thank D. Benedetti and J. Henson for sharing the Monte Carlo data underlying their work [37] with us. 
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„fit 
% 




(AP,)2 


70k 


0.7 X 10-^ 


7.5 X 10^5 


0.680 


100k 


8.8 X 10-5 


39.5 X 10-5 


0.318 


200k 


13 X 10^5 


61 X 10-5 


0.257 



Table 1: Initial conditions (7g*,AQ* for the RG-trajectory providing the best fit to the Monte 
Carlo data [17|. The fit-quality (AD^)^, given by the sum of the squared residues, improves 
systematically when increasing the number of simplices in the triangulation. 

(i) First, we numerically construct a RG-trajectory gkigo, Aq), ^kigo, ^o) depending on the initial 
conditions ^O; •^0; by solving the flow equations (|4.10p . 

(ii) Subsequently, we evaluate the resulting spectral dimension I?^^^(T; go, Aq) of the correspond- 
ing effective QEG space-time. This is done by first finding the return probability P{T; go, Aq), eq. 
(|3.15p . for the RG-trajectory under consideration and then substituting the resulting expression 
into ()2.5p . Besides on the length of the random walk, the spectral dimension constructed in this 
way also depends on the initial conditions of the RG-trajectory. 

(iii) Finally, we determine the RG-trajectory underlying the CDT-simulations by fitting the pa- 
rameters go, Aq to the Monte Carlo data. The corresponding best-fit values are obtained via an 
ordinary least-square fit, minimizing the squared Euclidean distance 

500 „ 
{AV^f^Y. , (5.1) 

T=20 

between the (continuous) function 2?P^*^(T; go, Xo) and the points 2?^^"'"(T). We thereby restrict 
ourselves to the random walks with discrete, integer length 20 < T < 500, which constitute the 
white part of Fig. |4]and correspond to the regime (iii) discussed above. 

The resuhing best-fit values go',\f for the triangulations with = 70.000, N = 100.000, 
and N = 200.000 simplices are collected in Table [TJ Notably, the sum over the squared residuals 
in the third column of the table improves systematically with an increasing number of simplices. 
By integrating the flow equation for g{k), X{k) for the best-fit initial conditions one furthermore 
observes that the points (7g*, Ag* are actually located on different RG-trajectories. Increasing the 
size of the simulation N leads to a mild, but systematic increase of the distance between the 
turning point T and the GFP of the corresponding best-fit trajectories. 

Fig. [5] then shows the direct comparison between the spectral dimensions obtained by the 
simulations (blue curves) and the best-fit QEG trajectories (green curves) for 70k, 100k and 200k 
in the upper left, upper right and lower left panel, respectively. This data is complemented by 
the relative error 

^ v9^''{T;gl\Xl^)-Vf^^{T) 
P?^«(r;5g^SAfi*) 

for the three fits in the lower right panel. The 70k data still shows a systematic deviation 
from the classical value T)s{T) = 3 for long random walks, which is not present in the QEG 
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Figure 5: Comparison between the 3-dimensional CDT data-sets 70k (upper left), 100k (up- 
per right), and 200k (lower left) obtained in [U] (blue curves) and the best fit values for 
pQEG^ji. ^fit^ ^fit-j (igpggjj curves). The relative errors for the fits to the CDT-datasets with 
N = 70.000 (circles), N = 100.000 (squares) and N = 200.000 (triangles) simplices are shown in 
the lower right. The residuals grow for very small and very large durations T of the random walk, 
consistent with discreteness effects at small distances and the compactness of the simulation for 
large values of T, respectively. The quality of the fit improves systematically for triangulations 
containing more simplices. For the = 200k data the relative error is ~ 1%. 



results. This mismatch decreases systematically for larger triangulations where the classical 
regime becomes more and more pronounced. Nevertheless and most remarkably we find that 
for the 200k-triangulation that e < 1%, throughout. All three sets of residues thereby show a 
systematic oscillatory structure. These originate from tiny oscillations in the CDT data which 
are not reproduced by 2??^'~^(T). Such oscillations commonly appear in systems with discrete 
symmetries [57] and are thus likely to be absent in the continuum computation. As a curiosity, we 
observe that the QEG result matching the most extensive simulation with N = 200k "overshoots" 
the classical value ^^(r) = 3, yielding 2??^^(r) > 3 for T > 450. At this stage, the RG- 
trajectory is evaluated outside the classical regime in a region of theory space where the Einstein- 
Hilbert approximation starts to become unreliable. It is then tempting to speculate that larger 
triangulations may also be sensitive to quantum gravity effects at distances beyond the classical 
regime. 

We conclude this section by extending v'^^'^ {T; g^^ , X^^) obtained from the 200k data to the 
region of very short random walks T < 20. The result is depicted in Fig. O which displays 
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Figure 6: Comparison between the spectral dimensions obtained from the dynamical triangulation 
with 200k simplices (blue curve) and the corresponding P?^^(T; <^q*, Aq*) predicted by QEG 
(green curve). In the latter case, the scaling regime corresponding to the NGFP is reached for 
log(T) < —40, which is well below the distance scales probed by the Monte Carlo simulation. 

Vf^^{T) (blue curve) and P?^^(r; gg*, A{j*) (green curve) as a function of log(r). Similarly to 
the four-dimensional case discussed in Fig. [3l the function P?^*^(r; ^q*, Ag*) obtained for d = 3 
develops three plateaus where the spectral dimension is approximately constant over a long T- 
interval. For successively decreasing duration of the random walks, these plateaus correspond to 
the classical regime P^^^(T) = 3, the semi-classical regime where 2?^^^(T) ~ 1 and the NGFP 
regime where pP^^(T) = 3/2. The figure illustrates that Pp^^(T) probes the classical regime 
and part of the first crossover towards the semi-classical regime only. This is in perfect agreement 
with the assertion [17] that the present simulations do not yet probe structures below the Planck 



6 Discussion and conclusions 

In this work we analyzed the fractal properties of the effective space-times arising within Quan- 
tum Einstein Gravity (QEG) formulated in the continuum by means of the gravitational average 
action. These effects are, to some extent, encapsulated in the spectral, walk, and Hausdorff di- 
mension seen by a fictitious diffusion process set up on the effective space-times. Most remarkably, 
these generalized dimensions are found to depend on the length of the random walk, indicating 
that the effective QEG space-times possess a multifractal structure. In particular, we established 
the possibility of a "low energy fractality" which occurs already well below the asymptotic scal- 
ing regime governed by the non-Gaussian (UV) fixed point (NGFP) and is thus unrelated to 
Asymptotic Safety. 



scale. 
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In Sections [3] and [4] we studied this multifractal structure within the flat-space Einstein- 
Hilbert approximation to QEG. Thereby it turned out that the effective QEG space-times are 
comparable to standard fractals in the sense that the relation ()3.3ip is satisfied on all scales. Their 
Hausdorff dimension is constant and equal to the topological dimension of the (background) space- 
time. Thus the fractal properties do not originate from the QEG space-times "loosing points" 
at short distances but rather represent a genuine dynamical effect of quantum field theory. In 
contrast to the Hausdorff dimension the spectral dimension and the walk dimension seen by the 
diffusion process depend on the diffusion time T. In Fig. [3] we identified three regimes in which 
these generalized dimensions are constant for a wide range of scales. These are connected by 
short crossovers. For long random walks, QEG space-times have the same spectral properties 
as classical flat space, i.e., the diffusion process is regular, = 2, and the spectral dimension 
matches the canonical dimension Vg = d. Moving towards shorter diffusion times one encounters 
the semi-classical scaling regime where = 2 + d, Vg = 2d/{2 + d). For infinitesimal random 
walks r — )• 0, the properties of the effective space-times are controlled by the NGFP and we 
obtain 2?^ = 4 and Vg = d/2. On both of the latter plateaus the random walk is recurrent. 

While the results concerning the NGFP regime are exact and follow directly from the very exis- 
tence of the fixed point, the "low energy" properties of the function T>s{T) rely on the applicability 
of the approximate effective field equation (j3.3p . It may be used if the Einstein-Hilbert approxima- 
tion is sufficiently precise, and if no matter fields are coupled to gravity whose energy-momentum 
tensor could possibly dominate over the cosmological constant term. In fact, our derivation of the 
running effective dimensions made essential use of the proportionality {g^u)k oc 1/Xk which does 
follow from the vacuum Einstein-equation (|3.3p . but not necessarily from a more complicated 
field equation G^^^{{g)k) = ->^k{gfiu)k - ^'T^GkT^^^{{g)k) where G^^^ is a (higher derivative, etc.) 
generalization of the Einstein tensor, and Tjl^,{{g)k) the energy- momentum tensor of the matter 
system in the effective geometry described by {{g^u)k)- The multifractal properties at length- 
scales below the asymptotic NGFP regime explored in this paper can occur only if the curvature 
of space-time is governed by precisely the scale-dependence of A^. If, on the other hand, \k{g^lu)k 
is negligible as compared to a /c-independent term in GkT!],,({q)u) the effective geometry has no 
s,gn,LL scale dependence and hence no fractal feaWsfFofthi .eason we expect thlt in ,ea. 
Nature the onset of the fractal behavior is typically shifted towards considerably higher energy 
scales than expected from the pure gravity case discussed at the end of Section HI We hope to 
come back to this point in a future publication. 

In Sect. [5] we performed a direct comparison between the spectral dimension of the three- 
dimensional effective QEG space-times with the one measured in Causal Dynamical Triangula- 
tions (CDT) [47J. Notably, the best-fit RG-trajectory reproduces the CDT data with approxi- 
mately 1% accuracy for the range of diffusion times where the simulation data is reliable. The 
comparison of P^^"'"(T) with I>P^'^(T) in Fig. [6] furthermore establishes that the present Monte 
Carlo simulations neither probe the semi-classical plateau nor the scaling regime of the NGFP. 
This confirms the cautious remark in ref. |47j that present day Monte Carlo simulations are un- 

®For the sake of the argument we assume that = Gp^ is the conventional Einstein tensor here. 
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able to probe physics well below the Planck length. This assessment also resolves the apparent 
contradiction between the extrapolation result limj-^o ^?^^(^) ~ 2 and the QEG prediction 
limj'_j.o I'?^'"' (T) = 3/2: The fit function employed in analyzing the Monte Carlo data can not 
be reliable extrapolated to T = and misses essential structures. 

The same conclusion also holds true in four dimensions. Comparing the profiles of I??^*^(r) 
shown in Fig. [3] with the fitting functions used in the CDT [35] or EDT [48] simulations shows 
that all the Monte Carlo data points obtained are positioned on the infrared side of the turn- 
ing point of the RG-trajectories underlying the QEG effective space-times. They neither probe 
the semi-classical plateau or the scaling regime of the NGFP. Performing the extrapolation of 
limr_).o P^^"'" (T) based on the leading corrections to the classical regime does not reliably iden- 
tify the signature of a non-Gaussian fixed point in 'Ds[T). Depending on where the data is cut 
off, one obtains different tangents to the first crossover, which lead to widely different extrapola- 
tions for the value dg = T)s{T)\t=o- We believe that this is actually at the heart of the apparent 
mismatch in the spectral dimension for infinitesimal random walks reported from the CDT and 
EDT computations. 

In order to test our conjecture that the non-classical P^-values found in the simulations 
are due to the "low energy fractality" predicted by QEG in the Einstein-Hilbert approximation 
one could perform the following numerical experiment. One couples gravity to a matter field 
whose parameters are averaged such that the resulting S.G.T,^. .s approximately ^-independeniS 
and much bigger than \k{g^iu)k- For this case we expect that the effective geometry has no 
significant /c-dependence, so that the fractal features disappear, and the spectral dimension equals 
the classical one. On the other hand, the dimensional reduction implied by the NGFP cannot 
be destroyed in this way. Assuming the fixed point is also present in the gravity-matter system, 
it enforces a scale-dependence upon S-kG^T^^ and the other terms in the effective field equation, 
which is precisely such that the solutions behave as {g^u)k oc implying the dimensional 

reduction T>s = d ^ d/2. 

We close our discussion with the following remark. In |64] the data set displayed in Fig. H] 
has been fitted to an anisotropic gravity model of Hofava-Lifshitz type. Comparing the quality 
of this fit with the residuals displayed in Fig. [5l we observe that both theories fit the CDT data 
with approximately equal quality. This is in particular remarkable, if one takes into account 
that approximating T>s{T) within the flat-space Einstein-Hilbert truncation gives rise to two 
fit parameters only, instead of the three parameters of the anisotropic model. Thus it may be 
premature to conclude that the spectral dimension obtained from the triangulations unequivocally 
identifles the underlying continuum theory as an anisotropic gravity model. In this light it seems 
mandatory to improve the simulation data in order to pinpoint the pertinent fractal structures 
of space-time with sufficient accuracy to reliably identify the underlying continuum theory. 

^Of course, for this purpose one could also give a non-NGFP scale dependence to T*^ and try to observe its 
impact on {g^,y)k- (Note that Gk will be approximately constant for k < mpi). 
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